Senescence Can Explain Microbial Persistence 

I. Klapper ' 6 , P. Gilbert , B.P. Ayati d , J. Dockery a < 5 , P.S. Stewart 6 ' 6 



"Department of Mathematical Sciences, Montana State University, Bozeman, MT 
59717 

^Center for Biofilm Engineering, Montana State University, Bozeman, MT 59717 
c School of Pharmacy and Pharmaceutical Sciences, University of Manchester, Manch- 
ester, UK 

^Department of Mathematics, Southern Methodist University, Dallas, TX 75205 
e Department of Chemical and Biological Engineering, Montana State University, 
Bozeman, MT 59717. 



Corresponding author: Isaac Klapper, Department of Mathematical Sciences, Mon- 
tana State University, Bozeman, MT 59717. Tel. (406)-994-5231. Fax. (406)-994- 
1789. email: klapper@math.montana.edu. 



1 



Abstract 



It has been known for many years that small fractions of persister cells 
resist killing in many bacterial colony-antimicrobial confrontations. These per- 
sisters are not believed to be mutants. Rather it has been hypothesized that 
they are phenotypic variants. Current models allow cells to switch in and out 
of the persister phenotype. Here we suggest a different explanation, namely 
senescence, for persister formation. Using a mathematical model including age 
structure, we show that senescence provides a natural explanation for persister- 
related phenomena including the observations that persister fraction depends 
on growth phase in batch culture and dilution rate in continuous culture. 

persisters | senescence | batch culture | continuous culture 



1 Introduction 

It has been observed (Balaban et al., 2004; Bigger, 1944; Gilbert et al., 1990; Green- 
wood and O'Grady, 1970; Keren et al., 2004; McDermott, 1958; Moyed and Bertrand, 
1983; Sufya et al., 2003; Haarison et al., 2005; WiufT et al., 2005), dating to Bigger 
(1944), that many antimicrobials while effective in reducing bacterial populations 
are unable to eliminate them entirely, even with prolonged exposure. The surviv- 
ing cells, called persisters, may be small in number - Bigger (1944) reported less 
than 100 persisters out of 2.5 • 10 7 cells of Staphylococcus pyogenes after exposure to 
penicillin in some cases for example - but nevertheless are subsequently able upon 
removal of the challenging agent to repopulate. See Lewis (2001) for a general dis- 
cussion. This phenomonon has recently gained increased attention in the context of 
biofilms (Spoering and Lewis, 2001) where the persisting cells have the added pro- 
tection of a polymeric matrix, making them particularly dangerous (see, e.g., models 
of Roberts and Stewart, 2004, 2005). The protection of microbial populations by 
persisters, whether formed by senescence or some other mechanism, is expected to be 
enhanced in biofilms because of the propensity of biofilms to harbor slow-growing or 
non-growing cells. 

Persisters have a number of interesting characteristics: 

• Upon reculturing, persister cells enable repopulation. See, for example, Balaban 
et al. (2004). 

• Persisters do not pass their tolerance to their progeny, and progeny do not 
inherit any greater tendency to be persisters. That is, persisters do not appear 
to be genetic variants. See for example Keren et al. (2003) and Balaban et al. 
(2004). 
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• It is observed that persister cells apparently grow slowly or not at all in the 
presence of antimicrobial. See for example Balaban et al. (2004). 

• Persister cells demonstrate resistance upon exposure to multiple anti-microbial 
agents. That is, persister cells with respect to one antimicrobial agent can also 
be resistant upon exposure to different antimicrobial agents. As an example see 
Sufya et al. (2003), where survivors (from an E. coli batch culture) of a tetra- 
cycline challenge were also tolerant to ciprofloxacin and quaternary ammonium 
compound. 

• Bacterial cultures demonstrate biphasic killing patterns in response to antimi- 
crobial challenge. It has been suggested that this plateauing is a consequence 
of the presence of persisters. See for example Balaban et al. (2004). 

• In continuous culture experiments, persister fractions are observed to increase 
with decreasing dilution rates, see Sufya et al. (2003). 

• All other things being equal, a culture with an initially richer growth medium 
will produce more persisters, see for example Bigger (1944). 

• One of the more puzzling observations is that changes in the population fraction 
of persister cells are growth-phase dependent - generally, though not always, 
population increases do not occur until the later stages of the log phase or 
even until the stationary phase. For example, Keren et al. (2004) observed 
an increasing ratio of persister to total cell numbers over time in E. coli, P. 
aeruginosa, and S. aureus batch cultures. In Balaban et al. (2004) the authors 
went as far as to posit that differences in onset of persister cell increases between 
different strains of E. coli imply the existence of more than one type of persister. 

It been suggested that that persistence is a phenotypic phenomenon, (e.g. Balaban et 
al. 2004; Cogan, 2006; Roberts and Stewart, 2005). These references propose models 
in which cells are able to switch in and out of a protected, slow- or non-growing persis- 
ter state with probabilities that are dependent possibly on environmental conditions. 
In this paper, based on observations of microbial senescence (Ackerman et al., 2003; 
Barker and Walmsley, 1999; Mortimer and Johnston, 1959; Stewart et al., 2005), we 
instead propose an alternative simple mechanism that can explain all of the above 
mentioned properties. 

We make only three essential assumptions: 

1. bacterial cells age, 

2. older cells are more resistant than younger cells to antimicrobial challenge, and 



3. growth manifests as production of new cells. That is, upon division one progenic 
cell inherits the effects of age while the other does not. 
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We thus regard older cells to be the persisters. In fact, one might well label senescent 
cells to be a separate persister phenotype (as opposed to a youthful phenotype) though 
we do not stress this interpretation. 

The second assumption could be a consequence for example of decreased metabolic 
activity but the particular mechanism does not really matter for our results. (The 
genetic bases for senescence and persistence are just beginning to emerge; see for 
example Nystrom, 2005, Vasquez-Laslop et al., 2006, and Spoering et al., 2006.) In 
reference to the third assumption, we can interpret cell division as a new, youthful 
cell being born from an old one. 

We make the following additional non-essential specific assumptions for definite- 
ness: 

4. production rate of new cells decreases with age but remains greater than zero, 

5. cell death occurs at a constant rate, and 

6. a given concentration of applied antimicrobial will kill cells of sufficiently young 
age but will not affect older cells. 

These extra assumptions matter in the details but do not affect the qualitative results 
that we report with the exception that a non-zero growth rate (in assumption 0]) is 
necessary in order to enable persister cells to repopulate after antimicrobial applica- 
tion. For consistency with the notion of decreasing activity, we suppose that older 
cells grow more slowly than younger cells although this assumption is again not really 
necessary here; rather the focus is on senescence as a mechanism for resistance. 



2 Persister Model 
2.1 Functional forms 

We define b(t,a)Aa to be the size of the bacterial population between ages a and 
a + Aa at time t, and c(t) to be the growth media concentration at time t. Following 
observations reported by Stewart et al. (2005), we suppose that cells senesce at a 
linear rate in time; other functional forms of age increasing senescence could be used 
instead. In particular for a parameter A that we will call the senescence time, we 
define u{t, a)Aa, the media usage at time t of cells between ages a and a + Aa, by 

1l(r M n] -j a[(l-a/X)+t]c(t), a<\, 
U[C[t) ' a) ~ \ a£c(t), a>\, 

where a is a first-order rate constant and £ is a base usage factor. 

Similarly, we define g(t, a)Aa, the new cells born at time t from cells between ages 
a and a + Aa, by 

g(c(t),a) = (3u(c(t),a) 
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where (3 is a yield coefficient. Additionally, by assumption the death rate fi(a) 
takes the form 

/i(a) = fi 

where //q is a death rate constant. Finally, functional forms for antimicrobial appli- 
cation will be made below. 

We stress that these choices are made for simplicity and for consistency with 
available data from the literature. The only essential condition we require is that the 
applied antimicrobial agent exhibit less potency as cells age (see below). For example, 
age- dependence in u and g is unnecessary. Conversely, fi could be made functionally 
dependent, for example, on a, c(t) or A if so desired. 



2.2 Age Structure 

A mathematical description of age structure was first introduced in Lotka (1907), 
McKendrick (1926), and many such representations have been used since. Define 
b(t, a) to be the bacterial population density of age a at time t, and c(t) to be the 
amount of available media at time t. The equation governing b is then 

db db 

-^(t,a) + —(t,a) = -fx b(t,a) (1) 

where the term on the right-hand side reflects cell death. The death coefficient /io 
can be expected, in general, to depend on c and a. We suppress such dependence 
here for simplicity as it does not affect our results in a qualitative way. 

Equation is valid for a > 0. To obtain an equation for b(t, a = 0), i.e., the 
new cells at time t, we observe that such cells are "born" at time t from the existing 
population b(t, a), a > 0. For example, the subpopulation of cells between ages a and 
a + Aa produces g{c(t), a)b(t, a)Aa new cells. Summing then over the entire existing 
population at time t, we obtain 

POO 

b(t,0)= g(c(t),a)b(t,a)da. (2) 
Jo 

A similar equation applies for the media concentration: 

dc f°° 

— (t) = — / u(t,a)b(t,a)da. (3) 
dt Jo 

Equations ((H)-(JHJ) are supplemented by initial conditions 6(0, a) = bo(a) for some 
supplied function bo(a), and c(0) = cq for some supplied constant cq. For discussion of 
mathematical issues involved in age-differentiated systems such as the one considered 
here, see e.g. Cushing (1998), Webb (1985). 
We can write u(t, a) = as(a)c(t) where 



s(a) 



l-a/A + £, a < A, 
f , a > A, 
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(s(a) is the senescence factor) and so define a senescence-weighted total bacteria 
population 

B{t) = I s(a)b(t,a)da. 
Jo 

Then CO)-© become 

9b db . . A . 

m + da- = a>0 ' (4) 

b(t,0) = af3cB, (5) 
dc „ 

- = -ocfl, (6) 

with initial conditions 6(0, a) = b {a), c(0) = c . In the computations to follow we 
use Co = 1 kg/m 3 and 



&o(a) 



(10 2 cfu)A _1 (l -a/2A), a<2A, 
cfu/hr, a > 2A. 



With this choice, the total initial population is 10 2 cfu distributed linearly in age over 
the age interval [0,2A]. 

2.3 Antimicrobial Application 

We include the effect of applied antimicrobial consistently with assumption |U1 a given 
antimicrobial concentration d applied at time t to the bacteria population results in 
killing at rate 7 of sufficiently young (and hence susceptible) cells and does not effect 
older (and hence tolerant) cells. We assume for definiteness and consistency with 
our senescence assumption that tolerance age increases linearly with antimicrobial 
concentration. Other choices of age dependence, as long as they are monotone in age, 
can be made. 

Then equation (J3j) is replaced by 

db db 

dt + da = + ^( rf ' a )) fo ' a > 0, (7) 

where the rate hk is defined to be 

(j \ j 1, a < Sd, , s 

^ a ) = { o, a >5d. (8) 

Here 5 is an adjustable tolerance coefficient. For a given antimicrobial concentration 
d, small 5 means that cells become resistant at a relatively young age, and large S 
means that cells become resistant at a relatively old age. For numerical reasons, we 
slightly smooth the discontinuity in /i^ in the results reported below. 

Based on (JHJ) then, for a given antimicrobial concentration d, cells of age Sd or 
greater are designated as persisters. 
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2.4 Chemostat Model 

In addition to batch culturing, we consider a chemostat system for which equa- 
tions (jH)-© become 

db db 

_ + _ = - {m + T)bt a>0 , (9) 

b(t,0) = a(3cB, (10) 

% = -acB + T(C -c), (11) 
at 

where r is the chemostat dilution rate and Cq is the reservoir substrate concentration. 



2.5 Parameters 

The model we describe contains seven parameters a, (3, 7, 5, A, /xo, and £ (plus two, 
r and Co for the chemostat), see Table 1. 



Table 1: Variables and Parameters 


B 


weighted population at time t 


cfu 


b 


population per age at age a and time t 


cfu/hr 


C 


reservoir substrate concentration 


kg/m 3 


c 


substrate concentration at time t 


kg/m 3 


d 


antimicrobial dosage 


kg/m 3 


s 


senescence factor 




a 


substrate usage rate 


hr _1 -cfu _1 




cell yield coefficient 


cfu -m 3 /kg 


7 


cell killing rate 


hr -1 


s 


tolerance coefficient 


hr-m 3 /kg 


X 


senescence time 


hr 


fJ'O 


nominal cell death rate 


hr -1 


T 


chemostat dilution rate 


hr -1 




minimum substrate usage factor 





In fact, the only essential constraints on these parameters necessary for the results 
we report are that 7 > /i , i.e., the antimicrobial killing rate is greater than the cell 
death rate, that 5 > 0, i.e., cells become more tolerant with age, and that £ > 
so that persister cells are capable of repopulation. We set the parameter values as 
follows: 

• Growth and substrate usage parameters: for definiteness, we set the length of 
the log phase to be approximately 10 hr and the bacteria doubling time to 
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be approximately 0.75 hr, resulting in approximately 13.3 doublings in the log 
phase. Given an initial value of 10 2 cfu, we thus obtain approximately 10 6 cfu 
at the end of the log phase. These constraints require cqo. — 2 -13 ' 3 = 1CT 6 
and c$af3 = 1.33 = (0.75) -1 (in units as in Table 1). We set Co = 1 kg/m 3 
as a reasonable value for initial media concentration, and then use a = 1CT 6 
hr _1 cfu _1 , (3 = 1.33-10 6 cfu-m 3 /kg. (Roughly speaking, a -1 fixes the population 
at the end of the log phase and (a/3) -1 determines the length of the log phase 
via the doubling time.) The minimum substrate usage parameter £ is presumed 
to be small compared to 1, and needs to be larger than zero in order for persister 
cells to repopulate. Otherwise its value is unimportant. We set £ = 10~ 3 . The 
value of the cell death rate is unimportant as well (and could even be set to 
zero). We use /io = 0.05 hr" 1 . 

• Senescence parameter: Supposing significant senescence after about 16 gener- 
ations (Stewart et al. (2005) reports 1-2% per generation) with a cell division 
time of 0.75 hr, we then obtain A = 12 hr. 

• Antimicrobial parameters (except 5): As a typical antimicrobial dosage, we use 
d = 0.01 kg/m 3 . As a typical antimicrobial killing rate, we use 7 = 10 hr -1 . 

• Chemostat parameters: we allow r, the dilution rate, to vary. The other pa- 
rameter, reservoir concentration Cq, does not affect our results so is arbitrary. 

The only important quantity without an estimate is the onset age 5d at which bacteria 
become resistant to the antimicrobial (see equation (jHJ)). We assume that the onset 
age be equal to A, the senescence time, and thus 5 = X/d = 1.5 ■ 10 3 hr-m 3 /kg. 
Increasing (decreasing) 5 has the effect of increasing (decreasing) the onset age. 

2.6 Numerical Methods 

Equation (JJJ) along with (J3J), (0) are solved numerically. We use a moving-grid 
Galerkin method in age with discontinuous piecewise linear functions post-processed 
to cubic splines for the approximation space in age (Ayati and Dupont, 2002). We 
use a step-doubling method in time, Ayati and Dupont (2005). This combination was 
illustrated in Section 5, Ayati and Dupont (2002), using the same code as used for 
this paper. 

3 Results 

In Figure [T] we apply antimicrobial to a stationary phase culture, Fig. da), and 
to a log phase culture, Fig. IHb), and then "reculture" afterwards by removing the 
antimicrobial and adding fresh medium. Note the biphasic survival in each case: there 
is an initial sharp die-off immediately after antimicrobial application followed by a 
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Figure 1: Exposure of a batch culture to antimicrobial, (a) antimicrobial is applied 
during stationary phase at t — 17 and removed at t — 27 at which time surviving 
cells are recultured at the t = nutrient level, (b) antimicrobial is applied during log 
phase at t — 5 and removed at t — 15 at which time surviving cells are recultured at 
the t = nutrient level. 

second phase of die-off of older (persister) cells due to natural causes. We also note 
that the number of persisters in this example does not grow substantially until the 
stationary phase - this is a consequence of the delay between birth of new cells in the 
log phase and aging of those cells into persisters. That delay time depends on 5, our 
only free coefficient. Small S results in short delay in the generation of persisters (i.e., 
growth in persister numbers even in the log phase) while large 5 results in long delay 
in the generation of persisters. The growth at early times in persister numbers seen 
in both Fig. ^a) and Fig. ^b) is solely a consequence of our initial age distribution 
in the innoculum. Note for example that persister numbers actually decline during 
the second log phase in both Fig. |TJa) and Fig. Hfb). 

In Figure El we present a computational version of the persister elimination exper- 
iment conducted in Keren et al. (2004), see Figure 4 of that reference, qualitatively 
matching the results reported there. These authors observed that persister numbers 
could be driven downwards by frequent reculturing. In addition to serving as a com- 
parison test of the persistence model to experiment, our (computer) experiment serves 
to further emphasize the point that the senescence mechanism presented here repro- 
duces observational evidence that persisters are not formed in early log phase, and 
that by suppressing production of senescent cells, it is possible to suppress persister 
frequency as in Kener et al. (2004). 

In fact, the persister elimination experiment illustrated in Fig. |21 might be consid- 
ered as an approximation to a chemostat with roughly 2 hour turnover. So we also 
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Figure 2: Bacteria batch culture, diluted down to 10 2 organisms every two hours. 
Solid curve is the total population and dashed curve is the persister population. See 
Keren et al. (2004), Fig. 4. 



consider persister numbers in the chemostat system ([§])-([TT]) at steady state (see Sufya 
et al., 2003). In particular, by setting db/dt = in ((0)) and suppressing t dependence, 
we obtain 

b(a) = 6(0)e- (w,+r)a . 
Thus the persister population fraction P is given by 

p = Im 6(a) da = e -(M0+T)M 
/ °° b(a) da 

see Figure El We again note the qualitative match to the experimentally reported 
results in Sufya et al. (2003). In particular, note the rapid transition from high 
to low persister fraction separating the two regimes r smaller than and larger than 
the persister age 8d. For small r, i.e., slow dilution, persister fraction tends to a 
constant controlled by /io, namely e~ tJl ° Sd ~ 0.55 for the parameter values used here. 
For large r, i.e., fast dilution, persister fraction tends to zero. We remark that 
this characterization of persister fraction with respect to dilution rate, see Fig. El is 
independent of the details of substrate usage. That is, it is independent of equations 
()1U|) and (|lip. ^ is a l so independent of the form of the senescence factor. We note 
though that at sufficiently high dilution rate, washout of cells exceeds the maximum 
specific growth rate and thus the biomass concentration, including persisters, becomes 
zero in steady state. 
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Figure 3: Fraction of persisters in a steady state chemostat as a function of chemostat 
dilution rate. 



4 Discussion 

The explanation of persistence as a symptom of senescence is an attractive one. 
Tolerance (or resistance) due to age, possibly because of reduced metabolic activity, 
explains in a simple manner the persister characteristics listed in the Introduction: 

• Offspring of old cells are of course young, and able to quickly repopulate. 

• Offspring of old cells are young and metabolically active cells, and hence do not 
inherit low activity tolerance, i.e., are not themselves persisters. 

• Old cells grow slowly. 

• Metabolic activity is independent of choice of antimicrobial agent so that metabol- 
ically inactive cells can demonstrate resistance to multiple agents. 

• In the presence of an antimicrobial agent, non-persisters are killed quickly fol- 
lowed by slow die-off of persisters, i.e., biphasic behavior occurs. 

• Decreasing the dilution rate in a continuous culture allows more cells to age 
longer before washout, hence increasing the persister fraction. 

• Increasing the amount of initial growth medium allows more new cells to be 
generated and hence, once these new cells have aged, more persisters. 
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• Increase in persister numbers can appear to depend on growth phase. In par- 
ticular new cells are, obviously, young cells and so log phase growth does not 
affect persister numbers until enough time has passed for those new cells to age 
sufficiently Hence persister numbers do not significantly increase until later in 
the log phase or even the stationary phase (and may actually decrease before). 

All of these observations are qualitative properties of our model depending on 
assumptions (JI}-© and, we believe, essentially independent of the particular choices 
made in the other assumptions. 

A number of authors have previously suggested persisters to be switching pheno- 
typic variants (e.g. Balaban et al., 2004, Cogan, 2006, Kussell et al., 2005, Roberts 
and Stewart, 2005, Sufya et al., 2005, Wiuff et al., 2005), that is, that persisters are 
cells with the same genome but with different sets of genetic expression as "normal" 
cells, and that a given cell can switch back and forth between the two states. The 
resulting phenotype switching model consists then of cells transiting between persis- 
ter and non-persister phenotypes. Available data does not distinguish between the 
alternative persister explanations of senescence and phenotypic switching. However, 
phenotype switching models require transitions between phenotypes to be growth 
phase dependent in order to explain phase dependence of persister numbers - per- 
sisters are more abundant in late stationary phase than in earlier stages. Hence, it 
becomes necessary to introduce stage- dependent transition rates, perhaps by making 
those rates nutrient dependent. Further, observations indicate that persistence in a 
given colony is antimicrobial dependent, that is, the number of surviving persister 
cells can depend on the choice of antimicrobial. Thus a phenotype switching model 
would seem to require a special parameter set (indicating, possibly, a special persis- 
ter) for each distinguished antimicrobial. On the other hand, the senescence model 
can handle this issue simply by allowing for different antimicrobials to be effective 
up until different levels of senescence, amounting in our presentation to making the 
parameter 5 be antimicrobial dependent. So, while we are unable yet to distinguish 
between phenotypic and senescent explanations of persisters based on available data, 
we can argue that the mechanism of senescence provides a simpler explanation at 
least with regards to the phenomenon discussed in this paragraph. 

We remark in conclusion that since (i) the segregation of chromosomes between 
mother and daughter cell is non-random and favors the parent with respect to the 
original rather than the copy strand, and (ii) mutations accumulate with successive 
replications of the chromosome then senescence/persistence/conservation of the older 
cells, particularly at times of antimicrobial stress, will ensure that following catas- 
trophic stress of a population, repopulation will be from older, demonstrably success- 
ful "archive" cells. This would avoid evolutionary dog-legs and possible extinction. 
Furthermore, consequent reliance upon the younger, mutation-prone cells to populate 
the community would increase the chance of adventitous genotypes being available 
in times of stress, whilst maintaining the "fall-back" position of reincarnation for the 
archived persisters. 
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